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We investigate the effects of two-body currents on elastic electron-deuteron scattering in an ex- 
actly Poincare invariant quantum mechanical model with a null-plane kinematic symmetry. While 
calculations using single-nucleon currents as input produce good qualitative agreement with exper- 
iment, the two-body current that we investigate produces a good quantitative agreement between 
theory and experiment for all three elastic scattering observables. 
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I. INTRODUCTION 

In the one-photon-exchange approximation elastic electron-deuteron scattering observables are functions of the 
current matrix elements. ((1, d)P' dl /^|/ M (0)|(1, d)Pd, fi), where |(1, d)Pd, /i) are deuteron four-momentum eigenstates 
and J M (x) is the deuteron current density operator. 

In a relativistic quantum theory the deuteron eigenstates transform irreducibly with respect to a dynamical unitary 
representation of the Poincare group. In this work they are eigenstates of a two-nucleon mass operator constructed 
by adding "realistic" nucleon-nucleon interactions, such as Argonne V18 or CD Bonn [HQ, to the square of the 
non-interacting invariant mass operator, M 2 = Mq + Amv nn . With this mass operator the S-matrix elements of the 
resulting two-nucleon model are consistent with the nucleon-nucleon scattering data used to determine the interaction. 

The representation of the dynamics can additionally be chosen so the interaction v nn is invariant with respect to a 
non-interacting, "kinematic" subgroup of the Poincare group. Different choices of kinematic subgroup 0] are related 
by unitary transformations that do not change the two-body S-matrix but the choice of representation of the 
dynamics effects the representation of the current operator. In this work the interaction is taken to be invariant with 
respect to a subgroup of the Poincare group that leaves a plane tangent to the light cone invariant. This is called 
the null-plane kinematic subgroup. It is the largest possible kinematic subgroup and it has well-known advantages 
for treating electron-scattering problems. 

The hadronic current density operator is conserved and covariant with respect to the dynamical representation of 
the Poincare group. It is necessarily the sum of one-body and two-body operators. Cluster properties imply that 
one-body current operators, which are determined by empirical nucleon form factors, are conserved and transform 
covariantly with respect to the non-interacting representation of the Poincare group. In a dynamical model the 
one-body parts of the current are only covariant with respect to the kinematic subgroup. 

Two-body parts of the current operator are needed to restore full covariance and may include additional contribu- 
tions of a dynamical origin. Deuteron matrix elements of an exactly covariant hadronic current are constrained by 
the Poincare transformation properties of the current operator and deuteron state vectors. These linear relations can 
be used to express all of the current matrix elements in terms of three linearly-independent current matrix elements, 
which are linearly related to invariant deuteron form factors. This current is also covariant with respect to the null- 
plane kinematic subgroup. The kinematic subgroup is relevant because the one and two-body parts of the current 
operator are separately covariant with respect to the kinematic subgroup. This more limited symmetry allows all 
current matrix elements to be expressed in terms of four linearly independent current matrix elements, using only 
null-plane kinematic transformations. The full covariance implies an additional dynamical constraint among these 
four kincmatically independent matrix elements. In a dynamical model with a null-plane kinematic symmetry the 
additional dynamical constraint is rotational covariance [f| [|[ Q [H which was appropriately referred to as the "angular 
condition" by Leutwyler and Stern. 

Any kinematically covaraint current operator can be used to define matrix elements of a fully covariant current by 
evaluating any three linearly independent current matrix elements using the kinematically covariant current operator. 
The remaining current matrix elements are generated by the constraints implied by covariance, current conservation, 
and parity. This construction ensures that the fully covariant current and kincmatically covariant current operators 
agree on the three independent matrix elements and all matrix elements related to them by kinematic Poincare trans- 
formations. The resulting covariant current depends on the choice of the three linearly independent matrix elements, 
which are chosen among the four kinematically independent matrix elements related by the angular condition. 

Null-plane "impulse calculations" compute the three independent current matrix elements using the single-nucleon 
current operators. When the null plane is oriented so the "+" component of the momentum transfer, Q + = Q° + Q 3 , is 
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zero the independent matrix elements can be chosen to be matrix elements of the + component, / + (0) = 7°(0) +/ 3 (0) 
of the current. When the independent matrix elements are taken to be matrix elements of the + component of the 
current the one-body matrix elements arc (1) invariant with respect to null-plane boosts (2) the single-nucleon form 
factors exactly factor out of the integral for the current matrix element and (3) the momentum transfer to the nucleons 
is identical to the momentum transfer to the deuteron. 

When used with realistic nucleon-nucleon interactions and empirical nucleoli form factors, null-plane impulse calcu- 
lations give a good qualitative description of all of the elastic electron-deuteron scattering obscrvablcs. The results are 
not sensitive to the choice of realistic interaction, empirical nucleon form factors, or the choice of independent matrix 
elements used to construct the fully covariant current. Quantitative differences between the results of null-plane 
impulse calculations and experiment, particularly in the observables A and T 20 cannot be explained by experimental 
uncertainties or uncertainties in theoretical input. 

In order to account for the difference between theory and experiment it is necessary to introduce an additional 
two-body current operator that has a non-vanishing contribution to the set of independent current matrix elements 
used to generate the fully covariant current. In this paper we investigate the consequences of including an additional 
two-body current that has an operator structure motivated by "pair-current" contributions in covariant formulations 
of elastic electron-deuteron scattering. In this work the two-body current is constructed to be covariant with respect 
to the null-plane kinematic subgroup. 

Null-plane quantum mechanical models use Poincare covariant two-component spinors rather than Lorentz covariant 
four-component spinors. Dirac spinors, which are spinor representations of Lorentz boosts, transform between the two- 
and four-component spinor representations. In the four-component representation, interactions that have a 75 vertex 
include contributions that couple to the photon through a causal propagator that do not appear in the null-plane 
impulse calculation. In the two-component spinor representation this contribution can be included as a two-body 
operator in the current. These considerations determine the structure of our model two-body current. The operator 
is constructed to satisfy exact null-plane kinematic covariancc. It is extended to a fully covariant current in the same 
way that the one-body current in the "null-plane impulse" calculations is extended to a fully covariant current. 

The resulting model current, when added to the current generated by the empirical nucleon form factors, leads to 
a good quantitative description of all of the elastic electron-deuteron scattering observables. Ambiguities related to 
the implementation of the angular condition are investigated. 

Elastic electron-deuteron scattering has been studied by many authors [100 0(1 00 00 
[00101 [ID] 21] [22] using a variety of different methods and assumptions. In all cases the problem is to con- 
struct matrix elements of a conserved covariant current operator for a range of space- like momentum transfers, Q. 
The calculations discussed in this work are based on quantum mechanical models where the Poincare group is an 
exact symmetry of the underlying theory and the interaction has a null-plane kinematic symmetry. Exactly Poincare 
invariant quantum mechanical models using other choices of kinematic subgroups have also been applied to electron- 
deuteron scattering. Different choices of kinematic subgroup have dynamical consequences for the impulse calculations. 
Another class of calculations are based on quasipotential treatments. These calculations extract the current matrix 
elements from covariant amplitudes without directly using the underlying quantum theory. This is done by general- 
izing Mandelstam's [23[ method to extract current matrix elements from Green functions using Bethe-Salpeter wave 
functions. "Relativistic impulse approximations" in this formalism depend on the specific quasipotential reduction, 
but they generally include a pair-current contribution that is related to the two-body current used in this paper. 
Detailed reviews can be found in 0] 0]. 

The Poincare invariant dynamical models in this paper are formulated in a representation with a null-plane kine- 
matic symmetry Q 0] 126] . Unlike models based on null-plane quantum field theory 0] 0] @] , where few-degree of 
freedom truncations break rotational covariance, the class of quantum mechanical models that we consider are exactly 
rotationally invariant. The null-plane treatment of electron scattering has the following advantages: 

1. For electron scattering the momentum transfer, Q, is spacelike, so the orientation of the null plane can be chosen 
so the +- component, Q + := Q a + Q 3 , of the momentum transfer is zero. When Q + = all current matrix 
elements can be constructed from three independent matrix elements of the + component, I + (0) := 7°(0)-|-/ 3 (0), 
of the current: 

((l,d),P', M '|/ + (0)|(l,d),P,^> P:=(P + ,Pi,P 2 ). (1.1) 

where we have labeled the deuteron four momenta by their null-plane components, P, and the spins are null 
plane spins (defined in section 2). 

2. There is a three-parameter subgroup of Lorentz boosts that leaves the null plane invariant. Because these 
boosts form a subgroup, there are no Wigner rotations associated with null-plane boosts. It follows that matrix 
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elements of I + (0) with Q + = and delta-function normalized initial and final states arc independent of frames 
related by null-plane boosts: 

((l,d),P"V|/ + (0)|(1, d),P",M> = ((l,d),P',f,'\I + (0)\(l,d),P,n) (1.2) 

where P'" = AP', P" = AP and A is any null-plane preserving boost. This means that matrix elements of / + (0) 
are equal to null-plane Breit-frame matrix elements with the same null-plane spin magnetic quantum numbers. 
The null-plane Breit-frame matrix elements are defined to have momentum transfer perpendicular to the unit 
vector, (z), that defines the orientation of the null plane. This is consistent with the requirement that Q + = 0. 

3. A consequence of property two is that the one-body current matrix elements of / + (0) with the physical momentum 
transfer exactly factor out of the corresponding nuclear current matrix elements: 

(*,P + Q|i+(0)|*,P) = 

(Pio + Q|J i f (0)|p; ) J dpi-- -dp„(#,P + Q|pi-- -Pi + Q---p„)(pi---pi---p„|*,P). (1.3) 

This follows directly from (|1.2j) because covariance with respect to null-plane boosts can be used to remove 
the dependence of the constituent current matrix elements on the internal momenta without changing the spin 
sums or normalization. A proof can be found in [26[. This means that the matrix elements of the one-body 
contributions to I + (0) can be expressed as sums of products of null-plane Breit-frame nucleon matrix elements 
multiplied by functions of Q 2 that only depend on the initial and final wave functions. 

With this formalism it is possible to satisfy exact rotational covariance and the impulse contribution to the inde- 
pendent current matrix elements involves only experimentally observable on-shell single- nucleon matrix elements. 

This paper is organized as follows. In section two we introduce our notation and define the null-plane kinematic 
subgroup. In section three we construct a dynamical representation of the Poincare group with a null-plane kinematic 
symmetry that provides a realistic description of the two-nucleon system. In section four we review the experimental 
observables for elastic electron-deuteron scattering and their relation to current matrix elements. In section five we 
discuss the nucleon currents that are used to construct the one-body part of our current operator. In section six 
we use the Wigner-Eckart theorem for the Poincare group to realize the dynamical constraints on the current, and 
we discuss some of the ambiguities in the implementation of the angular condition. In section seven we define our 
dynamical two-body current. In section eight we summarize our results. Appendix A summarizes how we relate our 
two- body current to the pion-exchange part of the model interaction, and appendix B contains additional details 
related to how this is done with the AV18 interaction. 



II. NULL-PLANE KINEMATICS 

In this section we discuss our notation and the null-plane kinematic subgroup introduced by Dirac 0]. The null 
plane is the three-dimensional hyperplanc of points tangent to the light cone satisfying the condition 

{x\x + :=t + x-e 3 = 0}. (2.1) 

The null-plane components, x = (x + ,xjj of the four vector x 11 are 

a: ± :=t±x-e3 x_l = (x • ei, x • 62). (2.2) 

Four vectors x^ can be represented by 2 x 2 Hcrmitian matrices in the null- plane components of x M : 

X=(l + X X *- ) = X*** *± = xi + iX2 x* = \Tr(^X) (2.3) 

where ct m = (I, a) and a are the two x two Pauli matrices. Since the determinant of X is minus the square of the 
proper time, — x 2 , if A is a complex 2x2 matrix with unit determinant then the transformation 

X -> X' = AATA f (2.4) 
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defines a real Lorentz transformation, 



A"„ = ^Tr^A^At). (2.5) 



We use the notation A to represent both the 2x2 and 4x4 representation of Lorentz transformations; the implied 
representation is easily determined. Points on the null-plane can be represented by triangular matrices of the form 

Poincarc transformations 

X -> X' = AAA f + A (2.7) 
that leave the null plane invariant (i.e. preserve the form of (|2.6[) ) have the form 

A =(?l/°c) < 2 ' 8 > 

where a ^ 0, (3, and a± are complex and a~ is real. These transformations form a seven-parameter subgroup of the 
Poincarc group that leaves the null plane invariant. This subgroup is called the kinematic subgroup of the null plane. 
This subgroup includes the three-parameter subgroup of translations on the null plane, a three-parameter subgroup 
of null-plane-preserving boosts as well as rotations about the 3-axis. 

The null plane-preserving-boost that transforms a rest 4-momentum (m, 0) to a value p is the matrix- valued function 
of v = p/m := {p + /m^p 1 /m,p 2 /m) given by 

Mp/m) : = A /W : = ( v j^t 1/^5+ ) 1 («) 

Since the reality of a in (|2.8p is preserved under matrix multiplication, the null-plane boosts, (|2.9[) . form a subgroup 
of the Lorentz group. 

The action of a null-plane boost on an arbitrary four momentum vector is determined by using (|2.9|) in (|2.4|) . The 
resulting transformation property of the + and _L components of the four momentum is 

p + ^p + ' = v + p + p± -> pi = p ± + v ± p + . (2.10) 

Since p~ does not appear in ()2.10j) the three components p := (p + ,p±) arc called a "null-plane vector". The — 
component can be calculated using the mass-shell condition 

p- = ^^. (2.11) 

p+ 

The null-plane spin of a particle of mass m is defined d so that (1) it agrees with ordinary canonical spin in 



the particle's rest frame and (2) is invariant with respect to null-plane boosts, (j2.9|) . 

The null-plane boost (|2.9p differs from the more standard rotationlcss boost, which has the 2x2 matrix form 



A c (k/m) = e^P a , (2.12) 

where p is the rapidity 



p = k|p| cosh(|p|) = y / k 2 +m 2 /m sinh(|p|) = k/m. (2.13) 

The null-plane representation of the single-nuclcon Hilbcrt space, 7ii, is the space of square intcgrablc functions of 
the null-plane vector components of the particle's four momentum and the 3-componcnt of its null-plane spin: 



V^PjaO = (0\ m )>P)/4 



/"DC /• J 

/ / rfpj. V |^(p,M)| 2 <oo. (2.14) 

Jo Jk 2 .._ ■ 



The Poincare group acts irreducibly on single-nucleon states 

((j,m)p,/i|£MA,A)|V) 
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/ dp+> dp' ± ]T ^if,',A A ' A ^,rn),p',M, 



(2.15) 



where the Poincare group Wigner-2? function in the null-plane irreducible basis [26l | is 

^;:;; ;p .,/A.,r : = <(j, m ),p,p\u(A 1 A)\( J , m ),p / 1 ^) = 



6(p-A(p')) x I^D^,[R fw (A,p')]e^ a , 



D 3 , [R] is the ordinary 5/7(2) Wigner D-function and 



R fw {A,p') :=Aj 1 (p/m)AA f (p'/m) 



(2.16) 



(2.17) 



is a null-plane Wigner rotation. It is the identity matrix when A is a null-plane boost. 

Our model Hilbert space for the two-nucleon system is the tensor product of two single- nucleoli Hilbert spaces. 
The kinematic representation of the Poincare group on the two-nucleon Hilbert space is the tensor product of two 
single-nucleon representations of the Poincare group: 



n = n 1 ®n 1 u (a,a) :=Ui(a,a)®Ui(a,a). 



(2.18) 



The tensor product of two one-body irreducible representations of the Poincare group is reducible. It can be decom- 
posed into an orthogonal linear superposition of irreducible representations using the Poincare group Clebsch-Gordan 
coefficients in the null-plane basis, 



*((ii, m i),Pi 5 Mi, (j2,m 2 ),p 2 ,fJ,2) = 

^2 J {(ju m i),f>uVr, (j2,m 2 ),P2,P2\{j,k)J,s,P,fi)dPk 2 dk 1 i'((j,rn)J,s,P,^), 

where the Poincare group Clebsch-Gordan coefficient [2(| in the null-plane basis is 

(C?i,mi),pi,Mi; (j2,m 2 ),p2,V.2\(j,k),l,s,P,{i) = 



(2.19) 



<5(P-pi-p 2 ) 



S(k - fc(pi,p 2 )) 



9(P,k) 



1/2 



k 2 <9(pi,p 2 ) 

]Tn m (k( Pll p 2 ))/y( 2 ^ 



(^'M'd \^2\ s ^s){l,m,s,fi s \j,p) 



with 



d(P,k) 


1/2 r 


9( P i,p 2 ) 


= V 



'(p++p+K(kV 2 (k) 



i? 



( Wl (k)+ W2 (k))p+ P +' 

f c (k/mi) := A^ 1 (k/mi)A c (k/mi) i? /c (-k/m 2 ) := Aj 1 (-k/m 2 )A c (-k/m 2 ) 



(2.20) 

(2.21) 
(2.22) 



and (u(k),k) = k(pi,p 2 ) = Aj 1 (P /M )pi. The rotation i?/ c (k/m), called a Melosh rotation, [3(|[26[ transforms the 
null-plane spins so they rotate under a single representation of SU(2) which allows them to be coupled using ordinary 
SU(2) Clebsch-Gordan coefficients. The absence of Wigner rotations in (|2.20p is a consequence of the null-plane 
boosts forming a subgroup. The quantum numbers I and s are degeneracy parameters that label different irreducible 
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representations with the same mass and spin that appear in the tensor product. In (|2.20[) the two-body invariant 
mass Mq, which has continuous spectrum, is replaced by k 2 which is related to M by 

M = 2 -A 2 + m 2 . (2.23) 
In this basis two-nucleon wave functions have the form 

<KV,k),i,sP, t i) = {(j,k),i,8p,M) / dp+ / dP ± / k 2 dk V V V |^((j,fc),z, s ,p, M )| 2 <00 

(2.24) 

and Uo(A,A) acts irreducibly on these states: 

((j,k),l,s,P,(x\U (A,A)\i>) = 



pOO n 3 

/ dP+> dP , x ^T^°^[A,Amk),l,a,P' f ^) 



W=—i 



I p+i 



E e iP '°y— ^[^(A^VMo)^^, fc),/, S ,P',/i'|^) (2.25) 

where P = AP'. This basis is used in the formulation of our dynamical model in the next section. 

In the rest of this paper we use the following notation for covariant, canonical, and null-plane basis vectors for 
mass-m spin-j irreducible representation spaces which are related by 



\(j,m),p,fx) = |(j, m ) )P , M )^^^M = J2 |(j,m),p, 2 /)^i?i M [P /c (k/m)]. (2.26) 



III. DYNAMICS 

Dynamical models of the two-nucleon system in Poincare invariant quantum mechanics arc defined by a dynamical 
unitary representation of the Poincare group acting on the two-nucleon Hilbert space. The mass Casimir operator for 
this representation can be defined by adding a realistic nucleon-nucleon interaction v nn [3l| [26j to the square of the 
non-interacting two-nucleon mass operator as follows: 

M 2 = M$+4mv nn , (3.1) 

where m is the nucleon mass and for a dynamical model with a null-plane kinematic symmetry and v nn has the form 

{(f,k')J\ s \P\ f i'\v nn \(j,k)J, s ,P^)^S(P'-P)6 rj S^^k\l\s'\\v^\\kJ, s ) (3.2) 

in the non- interacting irreducible basis (|2.20p . The interaction v nn is restricted so M 2 is a positive operator and 
designed so M 2 commutes with the non-interacting operators j 2 ,P and j z and is independent of the eigenvalues of 
P and j z . 

Simultaneous eigenstates of M 2 , j 2 ,P, j z in the non- interacting irreducible basis have the form 

((/, k'),l', s', P', M '|(j, A), P, M ) = S(P' - P)S j , j 6 ll ,^ Xtj (k, I, s) (3.3) 
where the wave function <p\ t j{k 1 1, s) is a solution of the eigenvalue equation 

1 \J+ S '\ poo 

(\ 2 -4k 2 -4m 2 )<j)x J (k,l,s) = Y J Yj I 4m(fc,;,s|KJ|fc',r,s')fc' 2 ( ifc> Aj (fc',r,s'). (3.4) 



s'=0l' = \j-s'\ 







The eigenf unctions, (f>xj(k,s), of this mass operator are also solutions to the non-rclativistic Schrodinger equation 
with interaction v nn , which can be seen by dividing both sides of (|3.4p by 4m. The deuteron is an even-parity bound 
state with j = s = 1 and parity limits the I' sum in (|3.4|) to I' G {0, 2}. 
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If {M 2 ,j, P,/x} have the same interpretation as {Mq,J, P,/x}, then it follows from (|3.2|) that the eigcnstates (|3.3[) 
transform irreducibly with respect to a dynamical representation of the Poincare group defined by replacing the 
eigenvalues of M in f- 0J -, ( [A, A] by the eigenvalues A of M: 

((j',kt),l',JP',i/\U{A,A)\(j,\),P,v l ) = 

J J2 (o^^),^^^p^M%^A),p^M , o^ , <(i.A),p / ^M // l^(A,i4)|(J^A),p,/i> = 

I*" =-3 

^(k^l^s'^f^J^A]. (3.5) 

Since the eigenstates (|3.3j) are complete, U(A,A), defined by (|3.5| . is a dynamical unitary representation of the 
Poincare group on the two-nucleon Hilbert space. When (A, A) is an element of the kinematic subgroup of the null 
plane, the Poincare Wigner 2?-function in the null-plane basis, , ^ [A, A], is independent of the mass eigen- 

value, A, and is thus identical to the Poincare Wigner 2?-function for the non-interacting irreducible representation. 
This shows that the null-plane kinematic subgroup defined by (|2.8|) is the kinematic subgroup of the dynamical 
representation (|3.5[) . 

The rotational covariancc in (|3.5[) is exact, but because T>^ ; ^ [A, A] depends explicitly on the mass eigenvalue 
A when A = R is a rotation, the rotations are dynamical transformations. 

The eigenfunctions of the interacting mass operator, M, are identical functions of k, I, s to the eigenfunctions of the 
non-relativistic Schrodingcr equation with interaction v nn . The identity 

Q±(H nr ,Ho nr ) = fl±(H r , Ho r ), (3.6) 

where fl±(H, Hq) are the scattering wave operators for the non-relativistic and relativistic system respectively, follows 
from the elementary calculation: 

n ± (H r ,H Qr )= lim e KPt+P-)t/2 e ^(P++P«)t/2 = Hm e iF-t' e -iP-f = 

t— >±oo t'— >±oo 

lim e ' t(M2+p o-L ) *' /p ^e" 4(M ° 2+p ^ )t ' /Po+ = lim e ^ 2 *" e -^o*" = 

t' — ► ioo t" — *±oo 

lim e i(^ 2 /rn+v nn )Arnt" +iArn 2 t" e ~i(\t 2 /rn)Amt'' -iArn^t" _ 
t" — '±00 

lim e i ( p2 /4m+k 2 /m+f nn )t'" e -i(P 2 /4m+k 2 /™)«'" 
' — ^±oo 

Q±(H nr , H nr) (3.7) 

where the limits are strong or absolute abelian [32j limits and the times in (|3.7j) were reparameterized as follows 

if" = Amt" = Amt' /P+ = 2mt/P+. (3.8) 

This proof uses kinematic symmetries that follow from (|3.2[) and the spectral condition, Pg > 0. It follows from eq. 
p.7p that relativistic and non-relativistic scattering matrices are identical functions of k: 

S(H nr , Ho nr ) = tt + (H nr , Ho nr )tt-(H nr , Ho nr ) = f2 + (H r , Ho r )n_ {H r , Ho r ) = S(H r ,Ho r ). (3-9) 

Since realistic nn interactions [H are constructed by correctly transforming experimental differential cross-section 
data to the center-of-momcntum frame and then fitting the solution of the scattering problem to the transformed 
data, these interactions can be used in (|3.I jl without modification. There is a small binding energy correction of 
about 1 part in 2000, which we ignore. The difference in the relativistic and non-relativistic wave function is due to 
the Poincare-group Clebsch-Gordan coefficients which are used to transform the wave function to the tensor product 
representation which is used in the computation of the current matrix elements. 

We use this method to construct the deuteron eigenstates and the representation of the Poincare group that we use 
to evaluate the deuteron current matrix elements. 
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IV. OBSERVABLES 

The differential cross section for elastic electron-deuteron scattering in the one-photon-exchange approximation, 



(27r) 4 e . .. . , , . , , . , , . (2tt) 3 

i, , 9 9 s | ( {3e , m e ) , p e , fx e 1 1 efl (0) [je , m e ) , \p e , /x e ) — — 
\J (p e ■ Pd) z ~ mjm^ q 



{{i,d),p' d ,^ d {0)\{i,d),p d ^ d )\H\ P ' d - Pd +p' e - Pe ) dp ' e d ff (4.1) 

is a quadratic function of the deuteron current matrix elements, ((1, d)p' d , /i d \I d (0)\(l, d),p d , fi d ). Covariance, current 
conservation, and discrete symmetries imply that only three of the deuteron current matrix elements are linearly 
independent [33|, which means that any deuteron clastic-scattering observable can be expressed in terms of three 
independent quantities for a given momentum transfer. 

Standard observables are the structure functions A(Q 2 ) and B{Q 2 ) and the tensor polarization T2q(Q 2 , 9) at 9i ab = 
70°. The quantities A(Q 2 ), B(Q 2 ) are determined from the unpolarized laboratory frame differential cross section 
using a Rosenbluth [34| separation: 



while T2o(Q 2 ,0) is extracted from the difference in the cross sections for target deuterons having canonical spin 
polarizations [i d = 1 and [j, d = at a fixed laboratory scattering angle (35| : 

TMQ 2 ,0) = V2 ^ iQ2 £\ n2 ^f Q2 ' 0) . (4.3) 

These experimental observables are related to the deuteron form factors @ HU [13 [H G (Q 2 ), Gi(Q 2 ), and 
G 2 (Q 2 ) by: 

A(Q 2 ) = G 2 (Q 2 ) + ^G 2 (Q 2 ) + G 2 (Q 2 ) (4.4) 

B(Q 2 ) = ^(l + r?)G 2 (Q 2 ) (4.5) 

T (C 2 6) = Gl[Q2) + ^ G °(Q 2 ) G2 (Q 2 ) + l^iW 2 )! 1 + 2 (! + g) tan 2 (^/2)] 

where 77 = j^p - and M d is the deuteron mass. 

The form factors are Poincare invariant quantities that are linearly related to current matrix elements. They are 
traditionally expressed in terms of canonical spin current matrix elements in the standard Breit frame, with the 
momentum transfer chosen parallel to the axis of spin quantization: 

G (Q 2 ) = d),Pl, 0|7°(0)|(1, d), P b , 0) + 2((1, d), Pi l|Z°(0)|(l,d), -P b , 1)) (4.7) 

Gi(Q 2 ) = -\[^((hd), Pi, WHom, d), P b , 0) (4.8) 

V v 

G 2 (Q 2 ) = ^((1, d), Pi 0|/°(0)| - (1, d), P b , 0) - ((1, d), Pi 1|/°(0)|(1, d), P b , 1)) (4.9) 

where P^ = (M d y/1 + rj, 0, 0, ^) an d P b = (M d y/1 + -q, 0, 0, — -y) and the normalization of G2{Q 2 ) follows the conven- 
tions used in 13811. 



9 



The O-momentum transfer limit of these form factors are related to the charge, magnetic moment, and quadrupolc 
moments 

lim G (Q 2 ) = 1 lim d(Q 2 ) = lim = -^=Q d . (4.10) 

Q 2 ^o Q 2 -o ' m n Q 2 ^o Q 2 2V3 

Current covariance, current conservation, and discrete symmetries can be used to express the matrix elements in 
(|4.7H4.9p in terms of matrix elements of the plus component of the current in the null-plane Breit frame, where the 
momentum transfer is in x direction and the spins arc null-plane spins. 

The canonical spin basis vectors in (|4.7H4.9[) are related to the null-plane basis vectors defined in (|2.14p by (|2.26|) 
which involves a change in normalization and a Melosh rotation (|2 . 22[) of the spins. 

V. NUCLEON CURRENTS 

The impulse current is the sum of single-nucleon current operators for the proton and neutron. For spin 1/2 systems 
that are eigenstates of parity there are two linearly independent current matrix elements (33j. They are related to the 
Dirac nucleon form factors F\(Q 2 ) and ^(Q 2 ) by 

(tin,m n ),p\v'\l»(0)\(j n ,m n )p,v}=u(p')r>*u(p) (5.1) 

where 

T» = rF 1 (Q 2 )-^-a^Q,F 2 (Q 2 ) (5.2) 
2m 

and u(p) is a canonical spin Dirac spinor. 

For spacelike momentum transfers with Q + = these form factors can be expressed in terms of the independent 
null-plane matrix elements 

<^ + (0)|i) :=(0„,m„),p',i|/ + (0)|(j„,m„),p,i)^F 1 (Q 2 ) (5.3) 
<^ + (0)| - \) := <(j„,m n ),p', i|/ + (0)|0„,m„), P ,-l) = -VfF 2 (Q 2 ) (5.4) 

where t := 

The Dirac nucleon form factors Fi(Q 2 ) and F 2 (Q 2 ) are related to the Sachs form factors by [39|] 

Fi(Q2) .g e (Q 2 ) + ^(Q 2 ) i (55) 

1 + T 

2 = Gm(Q2) ~ Ge(Q2) . (5.6) 

1 + T 

We consider recent parameterizations due to Bijker and Iachello [io| . Bradford, Bodek, Budd, and Arrington [4l| : 
Budd, Bodek, and Arrington [i^]; Kelly [43| and Lomon [44] 4 5| l46ll . These parameterizations all determine the 



proton-electric form factors using the polarization experiments |47j . The input to our calculations is the isoscalar 
form factors, which are the sum of the proton and neutron form factors. 

VI. DEUTERON CURRENTS: DYNAMICAL CONSTRAINTS 

The hadronic current density operator, I M (a;), transforms as a four- vector density under the dynamical representa- 
tion (|3.5|) of the Poincare group, 

U{k,A)I^{x)U\K 1 A) = {hT l y v I v (Kx + a). (6.1) 
The current operator must also satisfy current conservation 

<7 Mi 4PV(0)]_ =0, (6.2) 
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in addition to symmetries with respect to space reflections and time reversal. Here g^ v is the Minkowski metric with 
signature ( — V H — h). Current covariance and current conservation are dynamical constraints on the current operator. 
Using the identity 

((l,d),P',u'\I^x)\(l,d),P,p) = 

((1, d), P', u'\U^A, A)U(A, A)I»{x)U\k, A)U(A, A)\(l, d),P, v) (6.3) 

with current covariance (|6.ip and the transformation properties f|3 . 5[) of the deuteron cigenstatcs gives linear equations 
for each value of A and A that relate the different matrix elements: Current conservation leads to the additional linear 
constraints on the current matrix elements: 

9a ,(P' a ~ P a )<(l, d), P', ^(0)1(1, d),P,u) =0. (6.4) 

It is well-known that these constraints, when combined with space-reflection and time-reversal symmetries, [33| can 
be used to express any deuteron current matrix element as a linear combination of three independent current matrix 
elements. This is the Wigncr-Eckart theorem for the Poincare group applied to a conserved current. The invariant 
form factors or independent current matrix elements play the role of invariant matrix elements in the Wigner-Eckart 
theorem. 

If the coordinate axes are chosen so the + component of the momentum transfer is zero, which can always be 
done for electron scattering, then it is possible to choose the independent matrix elements as matrix elements of the 
"+" component of the current density, I + (0) 0. Matrix elements of I + (0) in an irreducible null-plane basis with 
normalization ()2.14j) are invariant with respect to the group of null-plane boosts. 

It follows from this invariance that the matrix elements of / + (0) only depend on the square of the momentum 
transfer and the individual null-plane spin components. There are 3x3 = 9 combinations of initial and final spins. 
Covariance with respect to the null-plane kinematic subgroup relates all of them to four kinematically independent 
current matrix elements. The one and two-body parts of the current operator transform independently under the 
action of the kinematic subgroup. The kinematic subgroup does not include the full rotation group. Rotational 
covariance [f| @ gives one additional constraint among the kinematically-covariant matrix elements, called the 
"angular condition" , reducing the number of independent matrix to three. It relates matrix elements of the one and 
two-body parts of the current operator. 

Because the constraints (|6.1j) and (|6.2|) are dynamical, it is a non-trivial problem to construct fully-covariant 
current operators. It is much easier to construct a set of deuteron matrix elements of a conserved covariant current 
operator. This can be done by first computing a set of three independent current matrix elements. The remaining 
matrix elements are generated from these three matrix elements using the constraints implied by (|6.1[) and (|6.4p . It is 
possible to do better by computing the three independent matrix elements using a model current that is covariant with 
respect to the null-plane kinematic subgroup. Operators that are covariant with respect to the kinematic subgroup 
are easily constructed. In this case the kinematically-covariant current and fully-covariant current will agree not only 
on the three independent currrent matrix elements, but also on any matrix elements generated from the independent 
matrix elements by the null-plane kinematic subgroup. The angular condition is needed to obtain all of the current 
matrix elements. It is a dynamical constraint that generates the fourth kinematically independent matrix element 
from the three independent matrix elements in a manner consistent with rotational covariance. The remaining current 
matrix elements can be generated from this fourth matrix clement using only kinematic covariance. This gives all 
deuteron matrix elements of a conserved covariant current operator. 

While it is straight forward to construct kinematically covariant current operators, there are many fully covariant 
current operators that agree with a given kinematically covariant current operator on different subsets of matrix 
elements. This leads to ambiguities when a kinematically covariant current is used to generate matrix elements 
of a fully covariant current. Model assumption are needed to eliminate these ambiguities. Specifically, a fully 
covariant current operator can be chosen to agree with the kinematically covaraint current operator on any set 
of three independent current matrix elements and all matrix elements generated from them using kinematic Poincare 
transformations. Because kinematic covariance cannot generate all of the current matrix elements, the resulting 
covariant current will depend on the choice of independent current matrix elements used to generate the full current, 
which is the model assumption. 

"Null-plane impulse approximations" choose the kinematically-covariant current to be the sum of the single- nucleon 
currents. For the deuteron the elastic scattering observables computed in the null-plane impulse approximation are 
not very sensitive to the choice of independent current matrix elements used to generate the full current. 

In this paper we use the same method to compute matrix elements of a dynamical two-body current operator. 
A kinematically covariant model two-body current is used to compute independent current matrix elements. The 
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remaining current matrix elements are generated using (|6.ip and (|6.4[) . As in the case of the null-plane impulse 
approximation, a complete specification of a fully-covariant model current requires both a kincmatically covariant 
current and a choice of kincmatically independent current matrix elements. 

If the deuteron eigenstates |(l,e£),P, fi) are given a delta- function normalization, then the matrix elements 

I+, y :=((l,d),P'^\I+(0)\(l,d),P,v) (6.5) 

with Q + = P +l — P + = and null-plane spins are invariant with respect to null-plane boosts. The four kincmati- 
cally independent matrix elements of / + (0) can be taken as I* 1 , 7^ , Iq and They are related by rotational 
covariance. The result is that only three linear combinations of these matrix elements can be taken as independent. 

The three independent linear combinations of the current matrix elements used to generate the fully covariant 
current in ref. Q were: 

^i+'oViVi + -i- ( 6 - 6 ) 

We refer to this choice of independent current matrix elements as choice I. These independent matrix elements can 
be distinguished by the number of spin flips between the initial and final states. 

Rotational covariance, or equivalently the angular condition, relates the difference Iyy — Iq to the matrix elements 
MM: 



its - ft,,, = ~YTv [v(I ^ + ~ + (6J) 

where 

7j := Q 2 /4M d 2 . (6.8) 

This relation is dynamical because r] involves the deuteron mass eigenvalue, Md- 

Direct computation of the difference I^y — Iq using the kinematically covariant current operator and the difference 
determined by the constraint (|6.7[) of current covariance 

A- := III - Ato + Wt,i + O - 2V2^i+o + lt,-x\ (6-9) 

gives a measure of the size of the dynamical contribution to the current operator that is generated by the rotational 
covariance constraint. If this difference is small compared to the size of independent matrix elements, then there will 
not be too much sensitivity to the choice of independent matrix elements; on the other hand, if this difference is large, 
there will be an increased sensitivity to the choice of independent current matrix elements. 

Both impulse or impulse plus an exchange current can be used to compute the independent current matrix elements. 
Our calculations, shown in fig. 1, indicate that A_ is larger when our model exchange current contributions are added 
to the independent matrix elements. This means that the presence of the exchange current enhances the sensitivity 
of the results to the choice of independent matrix elements. Because of this increased sensitivity, we investigate 
the impact of using different choices of independent linear combinations of current matrix elements to generate fully 
covariant current matrix elements on the sensitivity of the elastic electron-deuteron scattering observables. 

To motivate the other choices that we use note that any matrix element of a conserved covariant current in a set 
of deuteron eigenstates can be expressed in terms of a rank 3 covariant current tensor [36[ . The deuteron null-plane 
spin operator can be expressed in terms of the Poincare group generators as 

(OJ/) := -^Af 1 (P/M d )'> v e* a f h P a J*i (6.10) 

where P a is the deuteron four momentum operator, J' 37 is the deuteron angular momentum tensor, and Aj x {P/Md) v 
is a 4 x 4 Lorentz-transform- valued matrix of the operator P^ [2{| [26[ . The quantity (0, j/) is not a four vector because 
Aj 1 (P /MdY v is a matrix of operators. If the operator Af{P/Md)^u is applied to both sides of (|6. 10|) the result is a 
four vector, which up to the factor Md in the denominator, is the Pauli-Lubanski vector for the deuteron. 

It follows from (|6.10p that it is possible to replace the deuteron spins by four vectors by multiplying the current 
matrix element by the inverse of the boost used in (|6. 1 0[> ; the resulting vector is orthogonal to the four momentum. If 
this boost is applied to the initial and final deuteron states in a current matrix element then the result is a rank-three 
tensor density, T^(P',P) satisfying [2f3][H 

P p 'T^(P', P) = T^ a (P', P)P° = (P; - P^T^P', P) = 0. (6.11) 
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Formally this tensor density is related to the null-plane current matrix elements by 

T£(P', P) := A f (P'/M d Y P 0: H ^(d, P', v'\nO)\d, P, u)^O uj A f (P/M d y 



(6.12) 



1 

















1 


1 


i 





v/2 


V2 



where U j is the unitary matrix 

O vj = | 1 ] (6.13) 

that converts Cartesian components of vectors to spherical components and the i, j sums go from 1 to 3. The factors 
\J P + give |P, v)"/P+ a covariant normalization. 

The general form of this Lorentz covariant tensor density T^ a (P' , P) can be parameterized by the available four 
vectors and three invariant form factors Fi(Q 2 ), ^(Q 2 ) and ,F 3 (Q 2 ): 

T^(P',P) = l -(P' + PTE(P') pa [g^F l (Q 2 ) + T^-F 2 {Q^)]E{P)(s a 

- l -E{P') pa [Q a g^ - g» a QP]F 3 {Q 2 )E{P)p a (6.14) 

where 

E(P) aP =g aP -^ (6.15) 

is the covariant projector on the subspace orthogonal to P. The projectors are not needed when this tensor is 
contracted with vectors orthogonal to the final or initial four-momentum. 

We use this covariant representation of the current matrix elements to identify preferred choices of independent 
current matrix elements. Independent matrix elements can be constructed by contracting this current tensor with 
different sets of "polarization 4- vectors". These are 4-momentum dependent vectors v v (P) that are orthogonal P M . 
The contraction has the form 

v' af (P'YT^(P', P)v u {Py = \{P' + PYKf(P') ■ vuiP'W.iQ 2 ) + (<j(P ' ) 'y'" 6i(P)) F2(Q 2 )] 

- \[{<f{P') ■ Q)vu{PT - v af (P'r(Q ■ v bl (P))]F 3 (Q 2 ). (6.16) 

In this case the three invariant scalar products 

v' af {P') ■ v bi {P) v' af (P')-Q Q-v bl (P) (6.17) 

can be used in (|6 . 1 7|) and (|6.12j) to relate the form factors to the current matrix elements. Three pairs of polarization 
vectors give three equations that can be solved to express Fi(Q 2 ), ^(Q 2 ) and F 3 (Q 2 ) in terms of the current matrix 
elements. Given the form factors the rest of the current matrix elements are determined by the current tensor. 

The form factors F 1 (Q 2 ), F 2 (Q 2 ), F 3 (Q 2 ) defined by (j6T2"j) and ((5Tlj) are related to the form factors Gi(Q 2 ), 
G2{Q 2 ), G 3 (Q 2 ) [38l | which are defined in terms of canonical-spin Breit-frame matrix elements with momentum transfer 
in the 3 direction: 

((1 + 27?)^ + 477(1 - v)F 2 + 277F3) = c ((1, d), ^, 0|J°(0)|(1, d), 0) c = G (Q 2 ) + V2G 2 (Q 2 ) (6.18) 
F.iQ 2 ) = e ((l,d), |, l|I°(0)|(l,d),-|, l) c = G (Q 2 ) - -^f 2 {Q 2 ) (6.19) 
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F 3 (Q 2 ) = y| <(l,rf),§,-l|/ 1 (0)|(l,d),-|,0) c = G 1 (g 2 ). (6.20) 

The enhanced sensitivity to the choice of linearly independent matrix elements when the exchange current is 
included, as shown in fig. 1, suggests that in addition to the choice (|6.6p of linearly independent matrix elements used 
in Q, other choices should be considered. 

Frankrurt, Frederico, and Strikman (48j defined independent current matrix elements using polarization vectors 
(|6.16[) constructed from the last three columns vi c (P b ), v 2c (P b ), V3c(P b ), of the canonical boost in the null-plane Breit 
frame 

/ -Vv o o\ 

(MAW(A)«3c(A)) = io- ( 6 - 21 ) 

V o 01/ 

These are automatically orthogonal to the momentum. Three independent current matrix elements are extracted from 
the contractions v lc (Pl)PT+ (i* P b )v lc {P b ) a , i; 2c (P b ')"T p + (P b ', P b )v 2c {P b y and «3c(i?)T+. (P 6 ', P> lc (P h ) CT . These 
can be computed consistently using a kincmatically covariant model current in all frames related to the null-plane 
Breit frame by null-plane kinematic transformations. Their preference for this choice is to avoid using the 3'3 matrix 
element, which is maximally suppressed in the infinite momentum frame. In [j| operators are classified as good, bad 
or terrible according to their scaling properties with respect to Lorentz boosts to the infinite momentum frame. The 
3'3 matrix element is terrible in this classification. We call this choice of independent current matrix elements choice 
II. 

We are indebted to Fritz Coester (49| for suggesting a third choice of independent current matrix elements that 
emphasizes the kinematic symmetries of the null plane. This choice uses columns of the null-plane boost as polarization 
vectors. 

/ P x /M d py/M d P°/M d -M d /P+ \ 
(vi f (P),V2f(P),v 3f (P))=i I J P pJ^ (6 . 22) 

\-P x /M d -py/Md P z /M d + M d /P+ I 

The independent combinations are v ls {P')PT+ a (P' , P)v lf {P) a , v 2 f(P') p T+(P' P)v 2f {P) a , and 
V3 C (P') p T^ a .(P' , P)vi c (P) a . When the first two terms are evaluated in the null-plane Breit frame the scalar 
products (|6.17| that relate the contractions to the form factors are mass independent, which is preserved under 
kinematic transformations. This means that the boost parameters that appear in the calculations of the one-body 
current matrix elements arc independent of the constituent masses. To construct a third independent linear 
combination that exploits this mass independence would require using matrix elements of P(0). Rather than using 
P(0), P3 is extracted using polarization vectors that are mutually orthogonal, with one of them also orthogonal to 
Q, which gives matrix elements that are independent of Pi and F 2 . This results in the same expression for P3 = Gi 
used in the other two schemes. We call this choice of independent current matrix elements choice III. 

These three choices lead to expressions for dcutcron form factors based on three different covariant current operators 
generated using three different choices of independent linear combinations of matrix elements of a kincmatically 
covariant model current / + (0). These different current operators lead to distinct expressions for the deuteron form 
factors in terms of the null-plane current matrix elements of I + (0). The results, expressed in terms of the standard 
deuteron form factors and null-plane Breit frame matrix elements of the + component of the current are: 
Choice I: 

(1 + ry)G (Q 2 ) = {\- |)(/ n + loo) + + (f - \)h-i (6.23) 



(l + r,)G 1 (Q 2 ) = / 11 +/ ( 



(10 



'1-1 



(1 - V) 
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(6.24) 



(1 + V)G 2 (Q 2 ) 



I 1 + 4 A 

too) H 7T~ -'10 



V2 



(2 + V)h 



(6.25) 



Choice II (Go(Q 2 ) of choice I replaced by): 
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(1 + V)G (Q 2 ) = (f + l)i£ a - ftfo + ^i£ + (^)At-i] (6-26) 
Choice III (Go(Q 2 ) and G2{Q 2 ) of choice I replaced by): 

Go(Q 2 ) = (1 + f )/£ + llt-i y Gi{Q2) (6 ' 27) 

G 2 (Q 2 ) = - - '/Gx(Q 2 )) (6.28) 



Since B(Q 2 ) only depends on G\{Q 2 ) which uses the same linear combination of current matrix elements for all three 
choices, the computation of B(Q 2 ), is unchanged. These choices lead to a B(Q 2 ) that is in good agreement with 
experimental data. 

The different linear combinations of null-plane current matrix elements used in choices I , II and III have a non-trivial 
effect on the scattering observables A(Q 2 ) and T2q(Q 2 , 9). Selecting one of these choices is a model assumption. 

VII. DEUTERON CURRENTS: DYNAMICAL EXCHANGE CURRENTS 

The realistic interactions [l[ 0] used to construct the mass operator (|3.1[) include a one-pion-exchangc contribution 
plus a short-range contribution that is designed so the two-body cross sections fit experiment. 

"Required" two-body currents, discussed in the previous section, are needed to satisfy the constraints of current 
covariance and current conservation. These constraints arc satisfied by directly computing a set of linearly- independent 
matrix elements and then using the constraints to determine the remaining current matrix elements. Even if the 
independent current matrix elements are computed using the one-body impulse current, the matrix elements generated 
by the constraints will have two-body contributions. The interaction dependence arises because the covariance and 
current conservation constraints involve the deuteron mass, which is the eigenvalue of the dynamical two-body mass 
operator (|3.4|) . 

In addition to the two-body currents directly generated by covariance and current conservation, realistic interactions 
contain terms that involve the exchange of charged mesons, leading to the exchange of the nucleon charges, and the 
possibility of two-body currents associated with these exchanges. These currents give additional contributions to the 
linearly independent matrix elements of I + (0) that are needed to compute electron-scattering observables in null-plane 
quantum mechanics. The form of the dynamical exchange current contribution to / + (0) is model dependent, but we 
assume that the most important contribution is motivated by one-pion exchange physics. 

To motivate the structure of our two-body current consider the pseudoscalar pion-nucleon vertex 

£{x) = -ig n : *(x)^(x)t ■ <t>(x) : . (7.1) 

When this vertex is coupled to a photon-nucleon vertex with a causal propagator, the propagator connecting these 
two vertices can be decomposed into a sum of two terms involving the covariant spinor projectors A+(p) = it(p)u(p) 
and A_(— p) = —v(—p)v(—p), where u(p) and v(p) are Dirac spinors. This results in a sum of two terms with the 
following spin structures 

w(p)rMp'Mp'b 5 «(p") (7-2) 
-s(p)r'M-p'M-p'hMp")- (7.3) 

In this form these operator structures are expressed as 2 x 2 matrices in the nucleon spins. In Poincare invariant 
quantum models the first term is included in the one-body part of the current. The second term, which leads to the 
"pair current" in covariant formulations, is not included in the one-body part of the current in Poincare invariant 
quantum models. It can be included as a contribution to the two-body part of the current. To do this it is useful to 
convert the u-spinors to u spinors using 



«(-pM-p) = 75/Mp)u(p)/?75 



(7.4) 
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which replaces the operator structure in the second term by 

u(p)T»fa 5 u(p')u(p')pu(p"). (7.5) 

The factors (3 appearing in this expression are spinor space-reflection operators. Since space reflections do not leave 
the null plane invariant, they become dynamical operators in null-plane quantum mechanics. 

The canonical Dirac spinors can be converted to null-plane spinors by multiplying the spins by Melosh rotations, 
u f{p) = ' u >(p)D 1 '' 2 [R c f(p)]. These rotations cancel in the combinations that appear in the projectors A±(p), so the 
projectors can be expressed directly in terms of null-plane Dirac spinors. 

To construct a kinematically covaraint current we first make the model assumption that f3 should be replaced by an 
invariant operator that agrees with (3 in the rest frame of the initial (or final dcutcron). It requires /? — > —P d - j/M d . 

Because the null-plane boosts do not have Wigner rotations, the null-plane Dirac spinors satisfy the following 
covariance condition with respect to null-plane boosts A: 

u f (p) = SiA-^UfiAp) u f (p) = u f (Ap)S(A). 

This leads to the identity 

777 777 

-.u f (p)T+(-P d ■ j/M d ) l5 u f (p')u f (p')(~P d ■ 1 /M d )u } {p") 



TTL ~ ~ ~ ~ Tfl 

- 7 ==u } {Ap)T+{-AP d ■ 1 /M d ) l5 u f {Ap')u f {Ap')(-AP d ■ j/M d )u f (Ap") . (7.6) 

which defines a null-plane-boost-invariant current kernel that agrees with (|7.5p in the rest frame of the initial dcutcron. 
For the adjoint we use the rest frame of the final deuteron. 

The operator (|7.5[) factors into a product of 2 x 2 spin matrices, u(p)T M (— P d ■ j/M d )j 5 u(p') and u(p')(—P d ■ 
"f/M d )u{p"). The term u(p')(—P d ■ / y/M d )u(p") replaces u(p')75u(p") in one vertex of the one-pion-exchange inter- 
action, while the factor u(jp)T fl (—P d ■ 7 /Al d )j^u(p') has the appearance of a modified one-body current. 

To be consistent with our input two-body interaction we modify the one-pion exchange contribution of the model 
interaction by replacing the part of the phenomenological interaction that comes from u(p / )75u(p") by a modified 
interaction that comes from u(p')(— AP d ■ j/Al d )u(p") = u(p')(3u(p") in the deuteron rest frame. We then apply this 
modified interaction to the deuteron wave function and use kinematic null-plane boosts to transform the result to an 
arbitrary frame. 

The resulting matrix elements of the kinematically covariant two-body current has the following form 

((l,d),P', i/\I^(0)\(l,d),P,u) := 

f 1 P ■ 7 

/ ((M),P ,v Kin,^),^,^, (j n ,m n ),p 2 ,u 2 )(- — )u nf {p\)T^ lb {^^)u nf {pl)dp' 1 dp'[dp' 2 



X 



(tin,m n ),plvl{j n ,m n ),p 2 ,v 2 \U(A f (P/M d ))^\{lJ),P ,v) + (1 2) + h.c. (7.7) 

where in this expression we replaced the Dirac spinors with a covariant normalization used in projection operator 
A±(p) by spinors with a non-covariant normalization: 

Unfip) ■= J^Ufip), (7.8) 

to simplify the notation. We also remark that the projectors can be expressed in terms of canonical or null-plane 
spinors with covariant normalization, u c (p)u c (p) = Uf(p)uf(p). Details of the construction are discussed in appendix 
A and B. 

The quantity 

((1, d), P, v\(j„,m n ), pi, v\, (j n , m n ),p 2 ,v 2 ) (7.9) 

is the deuteron eigenstate in the tensor product basis, Po = (M d , 0, 0) is the rest frame value of the deuteron null-plane 
momentum, U(Af(P/M d )) represents a kinematic null-plane boost to the Breit frame. The modified pion-cxchangc 
interaction has the form (see Appendix A): 

(P / > k / ,i/ lJ i4|«, r |P,k > *i,*a ) 
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where ^(k — k') is the coefficient function of the one-pion-exchange contribution to the operator 

in the phenomenological interaction. For the Argonne V18 interaction ^(k — k') is extracted from the one-pion- 
exchange contribution to the tensor x (t\ ■ i~ 2 ) and spin-spin x (ti ■ r 2 ) parts of the interaction using the method 
discussed in [5(| and [HTj] . The extraction is discussed in Appendix B. The resulting interaction, v^{k — k'), differs 
from l/(m 2 + (k — k') 2 ) by the effects of the short-distance cutoff in the AV18 interaction. 

This exchange current is constructed to transform covariantly with respect to the null-plane kinematic subgroup. 
We add this current to the impulse current when computing the three independent matrix elements. These matrix 
elements are then used to generate the remaining current matrix elements using the constraints of current covariancc 
and current conservation as discussed in the previous section. 

Each term in the expression (|7.7|) for the exchange current matrix element can be represented as the product of 
( — 2m) with the modified current u n f(p[)T' 1 j 5 (j^)u n f(p'{) evaluated between a deuteron eigenstate and a pseudo- 
state defined by applying the rotationally invariant modified interaction « x in (|7. 10|) to the rest deuteron state and 
kinematically boosting the result to the Breit frame. This defines a kinematically covariant exchange current that we 
use to compute the independent matrix elements of 7 + (0). 



VIII. RESULTS 



The input to our calculation of the elastic electron-deuteron scattering observables is (1) the choice of nucleon form 
factors [HI [III [IB] [HI [H[ i (2) the choice of nucleon-nucleon interaction [l[ 0] and (3) the choice of independent 
linear combinations of current matrix elements used to generate the full current, 0[II[ [HI [12]. 

The extraction of proton-electric form factors based on polarization measurements compared with measurements 
based on the Rosenbluth separation were found to be inconsistent [4?| ; these inconsistencies have been explained [53| 
by including two-photon-exchange corrections in the Rosenbluth separation. This has led to a modification of the 
phenomenological parametcrizations of the proton electric form factor. All of the nucleon form factors that we use 
are consistent with the extractions based on the polarization measurements. 

Neutron electric form factor data is only available for a limited range of momentum transfers. The high-momentum 
transfer behavior of the different parameterizations is a consequence of different theoretical assumptions. This leads 
to some variation of the parameterizations for momentum transfers above Q 2 ~ l(GeV) 2 . Different parameterizations 
of the neutron electric form factor are compared to a dipole form factor in fig 2. The form factors that we compare are 
recent parameterizations given by Lomon J4f|, Budd, Bodeck and Arrington (BBA) [42|, Bradford Budd, Bodeck and 
Arrington (BBBA) [4l[, Kelly [43], and Bijker and Iachello (BI) [i(|. All of these parameterizations agree for Q 2 < 1. 
The curves in figure 2 are given as ratios to dipole form factors, which emphasize differences in the parameterizations. 
The input to our calculations are the isoscalar linear combinations of the nucleon form factors, F\n{Q 2 ) and F2n{Q 2 )- 
These form factors arc plotted in figures 3 and 4 for the same parameterizations that are compared in figure 2. These 
plots show no significant variation among the various nucleon form factors. Our calculations of elastic-scattering 
observables are not very sensitive to these differences. 

As shown in section VII, our model exchange current breaks up into a product of an effective "one-body" current 
and an interaction. The interaction is a modification of the one-pion exchange contribution to the tensor part of 
the phenomenological interaction. Most of our calculations are based on the Argonne V18 interaction. We extract 
the one-pion exchange contribution to the V18 potential by first discarding the short-range parts of the interaction 
that contribute to the spin-spin and tensor forces, then we use a method developed by Riska [50| and Schiavilla, 
Pandharipandc, and Riska [5l[ to extract the one-pion exchange contribution to the tensor force from the remaining 
parts of the interactions. This procedure is discussed in Appendix B. The one-pion-exchange potential that we extract 
differs from m a + (k_k') 2 ^ e enCec ts of the short distance configuration-space cutoff that appears in the Argonne V18 
interaction. The Fourier transform of the extracted interaction is compared to the pion-exchange potential without 
the cutoff in figure 5. The dotted curve includes the cutoff parameters that are used in the Argonne V18 interaction. 
The most important differences are for momenta above 1-2 /to -1 . 

Figures 6,7,8 show the three deuteron form factors, Go, G\ and G 2 with and without the exchange current included 
for the two-nucleon form factors (BI and BBBA) that have the largest high-momentum transfer difference in the 
neutron electric form factor. The independent matrix elements (|6.6p are calculated using the one-body parts of 
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the current and with the exchange current added. The remaining current matrix elements are determined by the 
constraints of current conservation and current covariance. The figures show that the addition of the exchange 
current contributions lead to an enhancement of Go above the minimum at Q 2 = IGeV 2 . The minimum of G\ shifts 
to the right, and G2 is enhanced. Except for momentum transfers Q 2 between 6 and 7 these form factors are not very 
sensitive to the different assumptions made about the neutron electric form factor. 

Data for the observables A are labeled Stanford Mark III [Hi], CEA [jtfLOrsay [HE], SLAC E101 [HI, Saclay 
ALS [H3, DESY [H, Bonnjp, Mainz J60[, JLab Hall C [H, JLab Hall A [gland Monterey Data for B arc 
labeled: SLAC NPSA NE4 (64^ Martinj6ff, Bonn [HI, Saclay ALS [HE], Mainz]^, Stanford Mark III 0. Data for 
T 20 are labeled: Novosibirsk-85 [13 [6|7Novosibirsk-90 [gg, Bates-84 0, Bates-91 JzH and JLab Hall C (z^. 

Calculations of A, B, and T 2 q using the independent matrix elements (|6.6[) (choice I) and the five parameterizations 
of the nucleon form factors used in figures 2-4 are shown in figures 9-11. While the null-plane impulse calculations give 
a qualitative understanding of the data, it is clear from these calculations that the null-plane impulse approximation 
is inadequate. 

Figures 12,13,14 show the effects of including the phenomenological pion-exchange current defined in section VII. 
We see that A and B provide acceptable fits to the data when the pion exchange current is included. The results are 
insensitive to the assumptions used in parameterizing the high-momentum transfer behavior of the neutron electric 
form factor. T20 is closer to the data, but it is still below the most recent Jlab hall C data [72j between Q 2 of .5 and 
2GeV 2 . 

The calculation displayed in figures 12,13,14 are based on the choice I of independent current matrix elements given 
by (|6.6[) . The presence of the exchange current (|7.7[) increases the sensitivity to the choice of the independent current 
matrix elements. This is shown in figure 1. What is plotted is the difference 1^(0) — Iq (0) with and without the 
exchange current using a direct calculation of the difference or by generating the difference using current conservation 
and current covariance. The difference between the dashed curve and solid curve shows that the required two-body 
contributions to the current in this difference is small when the independent current matrix elements are computed in 
the impulse approximation. Comparing the dotted and dash-dot curve indicates that much larger required two-body 
contributions to the current are needed when the exchange current contributions arc included in all matrix elements. 
This suggests that there will be a non-trivial sensitivity to the choice of independent current matrix elements used to 
generate the fully covariant exchange current. 

To test this we examined the two other choices, II and III, of independent matrix elements discussed in section VI. 
These methods relate form factors to independent current matrix elements by contracting different sets of polarization 
vectors into the current tensor. In both approaches there are preferred polarization vectors; in one case the vectors are 
chosen to minimize the dependence on matrix elements that arc maximally suppressed in the infinite momentum frame 
(P + — > 0), while the other choice minimizes the mass dependence in the contractions used to define the independent 
current matrix elements. Both choices were discussed in section VI. 

Figures 15,16,17, and 18 show the deuteron elastic scattering observables A(Q 2 ) and T 2 o(Q 2 ) for both choices, II and 
III, of polarization vectors. In both cases G\ is computed using the same linear combination of current matrix elements 
used for choice I. Since B only depends on Gi, B is identical for all three choices. All three choices of independent 
matrix elements give different predictions for A and T20. For choices II and III there is a mild enhancement of A at 
higher momentum transfers compared to choice of independent matrix elements given by (|6.6p . There is also a larger 
effect on the tensor polarization that brings the curve to within the experimental error bars. 

The result is that both choice II and III of independent current matrix elements give consistent results for the elastic 
scattering observables and they both provide a good description of the existing data over a wide range of momentum 
transfers. It is clear that there is a non-trivial sensitivity to the choice of independent current matrix elements when 
these results are compared to the corresponding results based on choice I. 

Another potential source of sensitivity to the input is the choice of nucleon-nucleon interaction. Any phase equivalent 
change in the nucleon-nucleon interaction is automatically accompanied by a corresponding change in the current 
operator. For interactions with a long-range meson exchange tail one might expect that the same data could be 
understood by simply adjusting the cutoff parameter. For typical soft interactions that arc useful in low-energy 
problems, one expects that a more significant modification of the current would be necessary. 

In figs. 19,20 and 21 we compare calculations of A(Q 2 ), B(Q 2 ), and T2o(Q 2 ) using the CD Bonn wave functions 
with and without the exchange current. In these calculations the exchange current is still based on the Argonne V18 
cutoffs. The calculations show the generalized impulse calculations and calculations where the exchange current is 
added to the impulse current, without adjusting the cutoff parameters. The calculations clearly show that there is 
more sensitivity to the choice of nucleon-nucleon interactions than to the choice of nucleon form factor. 

Good consistency with all experimental observables is obtained using the V18 interaction with nucleon form factors 
[ill ] and the choice of independent current matrix elements suggested by Frankfrut, Frederico and Strickman or 
Coester, discussed in section VI. These calculations are shown in figs. 22,23 and 24. For these choices the model 
exchange current explains the difference between the generalized impulse approximation and the experimental data. 
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Figure 25 shows the ratio of the experimental values of A(Q 2 ) to the calculated values using exchange current II. 
The data is presented on a linear scale to better illustrate the comparison between theory and experiment. The solid 
line represents the calculation with exchange current II while the diamonds show the ratio of exchange current III to 
exchange current II. There are no significant differences between the two exchange currents for momentum transfers 
below 4 (GeV) 2 . While our exchange current, which used the one-pion exchange part of the Argonne V18 interaction, 
required the numerical computation of Fourier transforms, the resulting interaction contribution to the exchange 
current differed very little from a simple static momentum-space pion exchange interaction, fig 5, which indicates the 
simplicity of our exchange current. 

Our results indicate that this simple exchange current is sufficient to provide a good quantitative understanding of 
elastic electron-deuteron scattering for a wide range of momentum transfers. The sensitivities to both the choice of 
interaction and the choice of independent matrix elements are the largest uncertainties in the calculations, and these 
uncertainties are all larger than the experimental uncertainties. For momentum transfers where data is available there 
is very little sensitivity to the uncertainties in the neutron electric form factors. 

Finally we compare the magnetic and quadrupole moments of the deuteron with and without the exchange current 
and using the three different choices of independent current matrix elements. The results are shown in table 1. Our 
calculations do not exhibit any sensitivity to the choice of nuclcon form factor, which are sufficiently well constrained 
by experiment at low momentum transfers. The exchange current contributions affect the results for all of the 
moments. In both cases they are closer to experiment than the moments computed using the generalized impulse 
approximation. The magnetic moment is in good agreement (to within computational accuracy) with experiment, 
while the quadrupole moment differs from the experimental result by a few percent. 

The conclusion of our research is that a simple exchange current motivated by one-pion exchange and the freedom 
to define a conserved covariant current operator by choosing a preferred set of independent current matrix elements is 
sufficient to provide a good fit to all three elastic scattering observables using Poincare invariant quantum theory with 
a null-plane kinematic symmetry. The model exchange current has the spin structure of a "pair" current, designed 
with a null-plane kinematic symmetry. The quality of the nucleon form factors has progressed to the point that our 
results are insensitive to the choice of nucleon form factor. 

While is it straightforward to include low-order pion-exchange physics in a more general class of models, the strategy 
for making the best choice of independent current matrix elements in a general class of electron-nucleus reactions 
requires more investigation. The principles used to derive choices I, II and III of independent current matrix elements 
all can be generalized to treat initial and final states with different spins. Whether there is one consistent set of 
principles that works universally for all reactions is not yet known. 

A second observation is that our model, with one of the two preferred choices of independent current matrix 
elements, provides a better description of all three observables than methods based on truncations of null-plane field 
theory or instant form relativistic quantum mechanics. Our model has features of both - unlike the instant-form 
model our model has the full-null plane kinematic symmetry with all of the advantages discussed at the beginning of 
this paper. Unlike truncations of a null-plane field theory, which emphasize cluster properties at the expense of exact 
Poincare invariance, our model is exactly Poincare invariant. We expect that the Poincare invariance constraint to be 
more important for momenta near or slightly above the deuteron mass scale, since the deuteron mass scale is involved 
in the implementation of the symmetry. 

This research provides a useful first step in trying to devise a more systematic treatment of model exchange currents 
in Poincare invariant quantum mechanics with a null-plane kinematic symmetry. It leads to a simple current that 
provides a significant improvement in all three elastic scattering observables when compared with the correspond- 
ing impulse calculations, however additional research is still needed in order to determine if these methods can be 
successfully applied to larger class of reactions. 

This work was performed under the auspices of the U. S. Department of Energy, Office of Nuclear Physics, under 
contract No. DE-FG02-86ER40286. The authors would like to express our gratitude to Fritz Coester who made many 
material suggestions that significantly improved the quality of this work. 

APPENDIX A: CURRENT CONSTRUCTION 

The steps motivating the structure of the model exchange current (|7.7[) are summarized in this appendix. We 
consider a pseudoscalar pion-nucleon vertex: 

C(x) = -ig* : ^(x)7 5 *(x)r ■ 4>{x) : (Al) 
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where — 2m^- is the pseudo-scalar pion-nucleon coupling constant and m is the nucleon mass. In what follows we 
use canonical Dirac spinors with a non-covariant normalization 



and use the notation 



TTl 

Wnc(p) := i/-r-r" c (p) (A2) 
w(p) 



9^ T T 

' V(k '' k) := (2ttV ml + (k'-ky-iO+- (A3) 
In this notation the tree-level one-pion-cxchangc transition amplitude is given by the rotationally invariant kernel: 

(k',//i,// 2 K|k,/ii,/4 2 ) := 

Wnc(k')75'«r lc (k)w 7r (k', k)u„ c (-k')7 5 u„ c (-k) (A4) 

where we have assumed a plane wave normalization (k'|k) = (5(k' — k). 

The tree-level one-pion-exchange transition amplitude in the presence of an external electromagnetic field, using 
the vertex (|A1[) with 

A»{q) := 7^4 J ^Mv^y, (A5) 

includes four terms, one of which is 

(p' 1 ,M / l ! P2,i«2l^(0) ,1 |pl,Ml ! P2,/"2) = 



(pi)r AI M„ c (r)u„ c (r)7 5 u„ c (pi) 



E12 -w(r) -w(p' 2 ) + i0+ 

^7r(P2! P2)Wnc(P 2 )75"r lC (P2) + 



u„ c (p 1 )r A1 w„ c (-r)w nc (-r)75Wnc(Pi)- 



where 



is the initial energy, 



and 



E12 - w(p' 2 ) + w(r) - i0+ 

«7r(P2 ! P2)w„ c (p 2 )75W n c(P2) H (A6) 

#12 =cj( Pi )+w( P2 ) (A7) 

r = Pi+p 2 -p 2 . (A8) 



=7^i(Q 2 )-7^^Q^2(Q 2 ) (A9) 

is the nucleon impulse current. The + • • ■ represents the contribution from the other three terms related by cither 
Hermitian conjugation and/or exchanging the proton and neutron. 

To motivate the structure our model exchange current we first evaluate this covariant expression in the rest frame 
of the initial two-body system so the pion-exchange interaction appears as a rotationally invariant function of the 
relative momenta. This allows us to relate the interaction part of this kernel to the interaction kernel in (|A4[) . Next 
we treat the initial energy of the two-body system, i?i 2 , as a parameter that can be expressed in terms of the mass 
of the initial state (deuteron) and the kinematically conserved momenta. We use the assumed null-plane kinematic 
symmetry of the interaction to express the energy denominators in terms of masses and null-plane kinematic variables; 
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E\2 — > \{P + + M p+ >x ), -E120 — > ^(-P + + M °p+' L )■ These are all model assumptions. In the rest frame of the initial 
deuteron P~ = Md which gives 



1 Af? - Ml k' Ml - Ml (k' M? - Af n 2 k' ) 
£12 - w(r) - w pU -> Bi2 - E 12a - — , ° { ' -> — ^ 2-^- « — ^ A10 



and 



1 il/ 2 + P i 1 A/ 2 
E 12 ~ w{p' 2 ) + w(r) - E12 -> + — — ) -» -(M d + -£) = M d » 2m (All) 



2 V P+ ' 2' M d 

Using ([ATO]) and (jATT|) in ((AG]) gives 

(p / lJ /x / lJ -k' jM ' 2 |7^(0)' 1 |k im) -k )M 2) = 

^»c(Pi)r A1 H„ e (k / ) M2 _ M 2 _ ■ 0+ Mnc(k / )75M»c(k)x 

u 7r (k',k)u„ c (-k')75W„ c (-k)- 

Wn C (pi)r A1 v nc (-k')-^-v„ c (-k')75U nc (k)x 
2m 

u w (k', k)u nc (-k')75W«c(-k) + • • • (A12) 

The terms in this expression are easy to interpret. The first two lines represent the product of a one-body current 
matrix element, the propagator, M 2_tS_ . 0+ , associated with equation (|3.4|) and the rotationally invariant kernel of 

^ d 

the pion-exchange interaction (|A4[) . This term is already included in the one-body contribution to the current matrix 
clement. 

The last two lines have a similar form except one of the u nc (k) spinors in the current and the rotationally invariant 
interaction kernel are replaced by v nc (— k) spinors. In addition, the propagator term is replaced by the factor 
l/2m. We find it convenient for computational purposes to split the v spinor terms that give the Dirac projector 
A_(— k) = — v c (— Vl)v c (— k) and to use 75/3 to convert the v nc (— k) to it rac (k). The two factors of (3 break kinematic 
covariance. We restore manifest kinematic covaraince by using this expression to define the current matrix elements 
in the rest frame of the initial system, and we transform to the Breit frame by requiring kinematic covaraince. 

Using 

Vnc(-P)v nc (-P) = 75/3u«c(p)Mnc(p)/?75 (A13) 

the last two lines of (|A12[) become 

u„ c (pi)r' 1 7 5 /3M„ c (k') 7 ^-x 

S nc (k')/3u(k)^(k'), (k))u nc (-k')76« n c(-k) + • • • . (A14) 
This separates into the product of an effective one-body current, 

u nc {v'i)^lsPu nc (\<!) (A15) 
a factor l/2m, and a modified one-pion-exchange interaction 

fZ„ c (k')/3u„ c (kK(k'), (k))u nc (-k')75W„ c (-k). (A16) 



The interaction term is identical to the rotationally invariant kernel of the one-pion-exchange interaction (|A4[) with 
the replacement: u„ c (k')75U nc (k) — > w nc (k')/3u„ c (k). This replacement preserves the rotational invariance of this 
kernel in the rest frame, but it is not covariant with respect to null-plane boosts. 
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This splitting has the advantage that the computation of the exchange current matrix element has the same structure 
as the computation of an impulse current matrix element with the current replaced by (|A15p and the deuteron wave 
function replaced by ^v\ip) where v is the modified interaction (|Al6j) . The input to this calculation is defined 
in rest frame of the initial deuteron. We get the null-plane Breit frame result by requiring kinematic covariance. 

The effective one-body current (|A15[) can be extended to a kinematically covariant operator that agrees with (|A15|) 
in the deuteron rest frame by replacing the (3 in (|A15j) by 

0->-p. 1 /M d (A17) 

This leads to a kinematically covariant modified current kernel 

«nc(pi)r^75^7^Wnc(k'). (A18) 

Kinematic covariance is all that this needed for a consistent computation of the current matrix elements. 

The modified interaction (|A16[) is rotationally invariant so when it is applied to the rest deuteron eigenstate the 
resulting pesudo-state has the same spin as the deuteron. This can be consistently defined in any other kinemati- 
cally related frame using null-plane boosts. This, along with the replacement (|A17[) restores the manifest kinematic 
covariance. 

The last step is to replace the modified interaction (|A16[) by the corresponding modified onc-pion-exchange part of 
a realistic model interaction. In a typical realistic interaction the pscudoscalar pion-exchange interaction is obtained 
by replacing the spinor terms in the rotationally invariant kernel (|A4|) by 

w„ c (k')75U„ c (k) -> - • ^~ k ' } (A19) 
2m 

This expression is obtained by retaining the leading term in a k/m expansion of the spinor term, which is normally 
justified because the one-pion exchange interaction also includes a high-momentum or short-distance cutoff. The 
net effect is that the resulting interaction, when included in the full nucleon-nuclcon interaction, provides a good 
description of the two-nucleon bound state and scattering observables. 

Expanding u(k')/3u(k) to the same order in k/m gives 1. This suggests that the modified interaction (|A16[) can be 
modeled by replacing the one-pion exchange contribution to the realistic interaction 

o-r ■ (k- k'K(k - k> 2 • (k' - k)n ■ r 2 (A20) 

by 

2mv n (k-k')er 2 • (k' - k)n ■ t 2 . (A21) 



The one-pion exchange interaction (|A4[) only contributes to the part of the tensor force in the nucleon-nucleon 
interaction that multiplies the isospin exchange operator T\ • t 2 . The tensor interaction also has contributions from 
the vector exchanges which also contribute to the spin-spin interaction. Riska [50| and Schiavilla, Pandharipande, and 
Riska [5l| introduced a method for isolating the pion-exchange contribution to the tensor force of a phenomenological 
interaction using linear combinations of the tensor and spin-spin interactions. The resulting interaction is 

v ps (k - k') = i(2«t(k - k') - v ss (k - k')) (A22) 

where v t and v ss are tensor and spin-spin contributions to the charge exchange part of the Argonne V18 interaction. 
We extract this from the pion-exchange contribution to the Argonne V18 potential. This is compared to mi + ^_ k ?)S 
in figure 5. The difference between these curves is due to the short-distance cutoff used in the AV18 interaction. 

When this interaction is applied to the deuteron bound state vector the result is a spin 1 "pseudo wave function" . 
The resulting "pseudo vector" can be defined in the Breit frame by requiring that it transforms covariantly with respect 
to the null-plane kinematic subgroup. This kinematic covariance ensures that the current kernel is kinematically 
covariant provided the pseudo-current is modified following (|A17[) . The kinematic covariance of the current is needed 
for a consistent calculation matrix elements of / + (0). 

We can now write the form of the Breit frame matrix elements of this one-pion exchange contribution to the 
exchange current: 

((1, d),P', fJl '\I+ (0)1(1, d),P,M> = 
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J2 J ((l,d),P\^,k%P''',l',s',^')k' 2 dk'dP'''x 
{(f,k'), P'"> l',s', fJ,'"\{j p , m p ),p' p ,fjf p , (j n , m n ),p' n , fj! n )dp' p dp' n X 

(Pp, Vp, Pn> ^ex-eff^Wp' m p)Pp, Mp: (jn, fn)Pn, fi n }dppdp n X 

(Op, m p ), p p , /Lt p , (i„, m n ),p n ,fj, n \(j, k)P", I, s, fi")k 2 dkdP"{(j, k)P",l, s, fj!'\P, fx, x) (A23) 
where the terms in this expression are the deuteron wave function in the free-particle irreducible null-plane basis 

((l,d)P',/i',d|(l,fc')^'^',P"',M"') - 8(P'-P"')5 j , f „6 fl ,^cl>*,(k',l , ,s') (f = s' = 1), (A24) 
the Poincare group Clebsch-Gordan coefficients in the null-plane basis 

(P'" > k',j',l', S '^'"\Pp^p,Pn^'n) = 



j( ^_ &; _ ft) *^P^^ y|(i|i|m , iSim<)x 

(s, m s |i, fi p , -,fJ. n )Y* mi (k(p^p;)) x 

^,[A- 1 (k/m)A / (k/m)]^,[A- 1 (-k/m)A / (-k/m)], (A25) 

the proton effective current 

<Pj )! / 1 p i Pn>/ i ll4- e //( )IPp^P i P'"/ 1 '>) = ( A26 ) 

5(p„-p„)u„ / (p;,M;)r^( ( 7)7 5 (^ 7 -?7 1 7 1 )«n/(pp,Mp), (A27) 
where the null plane Dirac spinors are related to the canonical Dirac spinors by a Melosh rotation 

w /(Pp:Mp)m = Ucippi^v'D^llk^iPp/T^kft-pp/m)], (A28) 
another Poincare Clebsch-Gordan coefficient 

(p p ,fx p ,p n ,fX n \P",k,j,l,s,fi") = (A29) 



_, _,^S(k'-k(p' p ,p' n )) I d(P,k) 



s(p»> - P ; - k) k , Vto x (A30) 

i 1 / 2 rA-Vw^A ^/^iin 1 / 2 



^,[A7 1 (k/m)A c (k/m)]^/ 2 < [A7 1 (-k/m)A c (-k/m)]x 



(s, m s |^, /ip, ^,n' n ) (j, mi,s, m s )Y lmi (k(p p , p' n )), (A31) 

and the pseudo wave function 

((fc,j),Z, S ,P", M "|P^,x) = 
S(P"-P) J 2j(j, //'|Z, to;, s, m 8 ) (a, m 8 |-, /x p , -, n n ) ^ 3/2 Y i* mi (k)rfkx 

(k - k') ■ o-^^fk' - k)n • T2dk'Y l , m ,{k')(s',m' s , \^,v' p , -,fx' n )(j,n, \, s',m' s , I'M^ujQs') (A32) 

where v ps is given by (|Al9jl . The full exchange current is the sum of the above quantity and the three other terms 
related by taking Hcrmitian conjugates or exchanging the particle that couples to the photon. The computation of 
the pseudo wave function is discussed in Appendix B. 
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APPENDIX B: CALCULATION OF PSEUDO WAVE FUNCTION 



The calculation of the pseudo wave function (|A32[) requires the one-pion exchange part of the tensor interaction. 
Riska [5(| and Schiavilla, Pandharipande, and Riska [5l|outlined a method to extract the pseudoscalar contribution 
to the tensor force as a linear combination of the radial coefficient functions in the Argonne V18 interaction. The 
method is based on the observation that both pseudoscalar and vector meson exchange contribute to both the tensor 
and spin-spin interaction in the static limit (50| . 

The Fourier transform of the interaction has the structure 

u n „(«) = Sl t v t (n 2 ) + SI ss v ss (k 2 ) H (Bl) 

where k := k' — k. 

fit := (cri • er 2 K 2 - 3<Ti • k<t 2 • k)(ti • r 2 ) (B2) 



Sl s :=K a (*i ■ o-aX-T! -r 2 ) (B3) 

and we have explicitly exhibited the tensor and spin-spin contributions to the interaction. 

The coefficient functions in (|B1[) are Fourier transforms of the coefficient functions the corresponding terms in the 
Argonne V18 interaction, 

4-7T f°° 

v ss (k 2 ) = — drr 2 v° r (r)(j {Kr)-l), (B4) 
K Jo 



Vt(K 2 ) = 



4tt 



drr v 7 '{r)j2(nr), 



(B5) 



where the factor 1/k 2 is due to the difference in the conventions used to define tensor operator in momentum and 
configuration space. 

These coefficient functions have contributions from both pseudoscalar and vector meson exchange. To separate 
them, following [50| . define the tensor and spinor operators and note that 

(o-i • kct 2 • «)((n • t 3 )) = --(Sh - Slss) (B6) 
while vector meson exchange gives the combination [Hoj ] 



This implies that 



((<ti x k) ■ (o- 2 x k))(ti • t 2 ) = \[2S1 SS + Sl t }. (B7) 



v t Sl t + v ss Sl ss = v ps ^(Sl ss - Sl t ) + v v -(2Sl s + Sl t ) = (B8) 



-{v v -v ps )Sl t + -(v ps +2v v )Sl s . (B9) 

Using these expressions we can isolate the pseudoscalar and vector contribution to the interaction using the following 
linear combinations of the tensor and spin-spin interactions: 

v ps = v ss - 2v t (BIO) 



v v = v ss + v t . (Bll) 

For the Argonne V18 interaction, which is used in our calculations, there remains a small difference between the 
pseudoscalar interaction calculated using (|B10[) and the contribution that comes directly from terms that can be 
directly identified with the one-pion exchange contribution. We only retain the one-pion exchange contribution to 
these terms, this still includes the short-distance cutoff used in the Argonne V18 r-space potential. 
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The interactions v° T (r) and v tT (r) have the form 

P (7 mo 
9 I lm± 



„<"•(,.) = £ 1 (-^ ) m Y( i i (h r) + 2m ± Y( t i±,r)} +v c (r) (B12) 



t^(r) = £ \( ^ ) m T(/x 0j r) + 2m±T( M± , r) |> + „*(r). (B13) 



9 1 \77l± 

Here y (/i, r) and T(/Lt, r) are the Yukawa and tensor functions with the exponential cutoff of the Urbana and Argoimc 
V\4 models: 

Y(n,r) = ^(l-e- cr2 ), (B14) 

T(/z, r) = (1 + i- + -J-)F( M r)(l - e^), (B15) 

The terms u c (r) and u'(r) are short-range phenomcnological interactions that we discard for the computation of the 
exchange current. 

In order to carry out the angular part of the integral in (|A32|) . we expand the pseudoscalar potential v ps (k — k') in 
partial waves: 

Vps (\k - k'|) = ]T Vl (k 2 , k' 2 )y ; : n (k)y im (k'), (Bie) 

where 



(k 2 ,k' 2 ) = 2tt J Pi(u)v ps {\/'k 2 + k' 2 - 2kk"u)du. (B17) 



We expand the vector k and k' in spherical harmonics 



k = fc(^(y 1 _ 1 ( k ) - ni(k)),*Y y (y x _i(k) + Fll(k)), JyFlo(k)) (B18) 



to obtain the following expression for the pseudo wave function: 

(c) iTii II 



X f d Jk",n!;,n' n ) (B19) 

-i /> I CO I 

55 T^r^E ^(k 2 ,k" 2 )y^(k")^(k)y iw (k) 

^ ' J !=0,2^,=-ii'=0^J=-i' 

x [fc"( J^(Fi-i(k") - y n (k")),i J — (Fi-i(k') + y n (k)), \/— Yio(k)) 



3 v , U l M ' v V 3 

- fc(y^(y 1 _ 1 (k)-Fn(k)),yy(y 1 _ 1 (k) + y 11 (k)),y T y 10 (k))] 

x (l4i\v\V>n) ( S P; A»p ) S n ,/i„|s ) /X 8 )(s,/i a ,/, /i; | j, m)-^!, ( k )^ (fc), 

The angular integrals are evaluated using: 



| dky^ I '(fc)y / i I (fc) = ^ w> (B20) 

JdkYfl'faY^&Yfcft) (B21) 
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After carrying out the integral on k: the expression for the pseudo wave function becomes 

X f d Jk"^';^' n ) (B22) 



t^/^E E E E vi^ 2 x 2 ,)Y^")Mk) 

1 ' J 1=0,2 Hi=-l V H'i=~l' 

X (s p , fjfp, S n , H n \s, (is) (s, ll s , I, (J,) 

Phi 



fc v (2 l^-t'+i) 1} (1, 0l *' 0|r ' 0)( ~ (1 ' -1, z ' w>) + (1, 11 u 



+ kJ (2 ±tl^_i± 1} (i, o, i, o\i>, o) (-(i, - i, i, w ,) - 



/ 47T 

+ \/y(M»«) [ fc"l r io(k")*i'i<W w 



The sum over I' includes only a finite number of values. It can only be or 1 for I = and 1, 2 or 3 when 1 = 2. This 
limits the sum on V to the first four partial waves. V = 0, 1, 2, 3. The final expression of the pseudo wave function can 
be written in the form: 



V 



' (£">Mp./4J =^2 I U'( k ")flt'(k",fJ'p,^n): ( B23 ) 



where fai (h" , [i' n ) are angle-dependent coefficients and Iw(k") are the relevant scalar quantities which have the 
following form: 

I w {k") = [ k 2 dk Vl (k,k")uv(k) (B24) 



or 

•J'\ I U'3 



I w {k") J k /3 dkvi(k,k") Ul ,(k). (B25) 
The allowed combinations of I and V pairs requires the following integrals Iw(k"): 

I 00 (k") = I ' v (k,k")u (k)k 2 dk (B26) 



I w (k") = J v 1 {k,k")u {k)k 3 dk (B27) 

I X2 {k") = J v l {k,k")iL 2 {k)k i dk (B28) 

I 22 (k") = J v 2 (k,k")u 2 (k)k 2 dk (B29) 

I 32 (k")= [ v 3 (k,k")u 2 (k)k 3 dk (B30) 



2G 

The total exchange current contribution is 

(P' d ,tx' d ,d\I+(0)\P d ,v d ,d) =J2j dk'Ji (B31) 

x ( k '> Mp, Mn)[«/(Pp, M P )r^(g)7 5 (-^ • t/Mj)«/(Pp, Mp)]*^ ( k ; A'p, Mn) 

+ X^?^(k, n' p , Mn)[%(Pp ; ^)(--P ' 7/Md)7 5 r^(g)w/(pp, Mp )]*p d Aid (k, M „)] 
+ [p n] 
where is the jacobian 

Ji - l w^) 1 (B32) 
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IM IM+Exchange IM(WSS) IM+MEC(WSS) 

Q d 0.2698 0.2752 0.270 0.275 

0.8535 0.8596 0.847 0.871 



TABLE I: deuteron magnetic and quadrupole moments evaluated in the impulse approximation and including the exchange 
current contribution. The values are the same using all six different nucleon form factor parameterizations and three combina- 
tions of independent current matrix elements. Argonne V18 potential is used in the calculation. The values labeled with WSS 
are from [J. The experimental values are 0.2860 ± 0.0015/m 2 [fl and 0.857406 ± O.OOOOOI^at [zl- 




fig. i: (color online) /+ + (0) — /qo(0) calculated directly compared to /++(0) — I, 
calculated by rotational convariance with and without exchange currrent I 
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Neutron Electric Form Factors 

E n dipole 




fig. 2: (color online) Neutron electric form factor parameterizations 
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fig. 3: (color online) Isoscalar nucleon form factor Fin(Q 2 ) 




fig. 4: (color online) Isoscalar nucleon form factor F2n{Q 2 ) 
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One-Pion Exchange Interaction 

comparison of OPE part of AVI 8 to OPE potential 
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fig. 5: (color online) One pion-exchange interaction compared to one-pion exchange 
part of AV18 interaction 
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fig. 6: (color online) Deuteron form factor Gq(Q 2 ); different nucleon form factors with 
and without exchange current I 
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fig. 7: (color online) Deuteron form factor G\(Q 2 ); ; different nucleon form factors 
with and without exchange current I 
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fig. 8: (color online) Deuteron form factor G^Q 2 ); ; different nucleon form factors 
with and without exchange current I 
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A(Q ): Generalized Impulse Approximation 

with different nucleon form factors 
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fig. 9: (color online) A(Q 2 ): Generalized impulse approximation with different nu- 
cleon form factors 
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B(Q ): Generalized Impulse Approximation 

with different nucleon form factors 
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fig. io: (color online) B(Q 2 ): Generalized impulse approximation with different nu- 
cleon form factors 
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fig. ii: (color online)T2o(Q 2 , 70°): Generalized impulse approximation with different 
nucleon form factors 
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A(Q ): Impulse + Exchange Current I 



with different nucleon form factors 
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fig. 12: (color online) A(Q 2 ): Generalized impulse approximation plus exchange cur- 
rent I 



41 




fig. 13: (color online) B(Q 2 ): Generalized impulse approximation plus exchange cur- 
rent I with differnet nucleon form factors 
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fig. 14: (color online) 72o(Q 2 , 70°): Generalized impulse approximation plus exchange 
current I with different nucleon form factors 
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A(Q ): Impulse + Exchange Current II 



with different nucleon form factors 




fig. 15: (color online) A(Q 2 ): Generalized impulse approximation plus exchange cur- 
rent II with different nucleon form factors 
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T 20 (Q 2 ,70°): Impulse + Exchange Current II 

with different nucleon form factors 
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fig. 16: (color online) 72o(Q 2 , 70°): Generalized impulse approximation plus exchange 
current II with different nucleon form factors 
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A(Q ): Impulse + Exchange Current III 



with different nucleon form factors 
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fig. 17: (color online) A(Q 2 ): Generalized impulse approximation plus exchange cur- 
rent III with different nucleon form factors 
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T 20 (Q 2 ,70°): Impulse + Exchange Current III 

with different nucleon form factors 
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fig. 18: (color online) T2o(Q 2 , 70°): Generalized impulse approximation plus exchange 
current III 
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A(Q ) : CD Bonn Interaction 

with and without exchange currents II and III 



i — i — i — i — i — i — i — i — r 



i — i — i — r 









Stanford Mark III 





CEA 


A 




□ 


SLAC E101 




Saclay ALS 





DESY 







> 


Mainz 


□ 


JLab Hall C 


o 


JLab Hall A 


X 


Monterey 




IM CDB -I 




1M+Exehange III 




IM+ Ex change II 



_l I I L_ 



0.5 



1.5 
Q 2 [GeV 2 ] 



2.5 



fig. 19: (color online) A(Q 2 ): CD Bonn with and without exchange currents II and 
III 
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fig. 20: (color online) B(Q 2 ): CD Bonn with and without exchange currents II and 
III 
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fig. 2i: (color online) T2o(Q 2 ,70°): CD Bonn with and without exchange currents II 
and III 
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fig. 22: (color online) A(Q 2 ), AV18, BBBA form factors with and without exchange 
current II and III 
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B(Q ): BBBA Nucleon Form Factors 

impulse and exchange currents II and III 
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fig. 23: (color online) B(Q 2 ) , AV18, BBBA form factors with and without exchange 
current II and III 



■52 




fig. 24: (color online) T2o(Q 2 , 70°) , AV18, BBBA form factors with and without 
exchange current II and III 
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2 2 

Ratio of A(Q ) experimental to A(Q ) calculated 

BBBA nucleon form factors, exchange current II 
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fig. 25: (color online) Ratio of experimental values of A(Q 2 ) to calculated values of 
A(Q 2 ) with exchange current II 



